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Abstract 

The Hartle-Thorne metric defines a reliable spacetime for most astrophysical purposes, for 
instance for the simulation of slowly rotating stars. Solving the Einstein field equations, we added 
terms of second order in the quadrupole moment to its post-linear version in order to compare 
it with solutions found by Blanchet in the frame of the multi-polar post-Minkowskian framework. 
We first derived the extended Hartle-Thorne metric in harmonic coordinates and then showed 
agreement with the corresponding post-linear metric from Blanchet. 

We also found a coordinate transformation from the post-linear Erez-Rosen metric to our ex¬ 
tended Hartle-Thorne spacetime. It is well known that the Hartle-Thorne solution can be smoothly 
matched with an interior perfect fluid solution with physically appropriate properties. A compar¬ 
ison among these solutions provides a validation of them. It is clear that in order to represent 
realistic solutions of self-gravitating (axially symmetric) matter distributions of perfect fluid, the 
quadrupole moment has to be included as a physical parameter. 
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1 Introduction 

In 1968, Hartle and Thorne (HT) [16l [26] proposed an approximate solution to the Einstein field 
equations (EFE) intended to represent the gravity field of neutron stars with mass, rotation and 
quadrupole moment as parameters. Berti et al. [1] compared the HT metric with the Manko [201 f2T] 
(exact solutions) and Cook-Shapiro-Teukolsky metrics |6j (numerical solution), and showed that it 
is safe to use the HT metric, since it gives excellent results even for the innermost stable circular 
orbits with fast spin periods. Moreover, the exterior HT metric can be smoothly matched with a 
physically reasonable interior one. This provides realistic models of stars and for this reason, it is 
often used to validate exact and approximate solutions of EFE. 

Stationary exact solutions of EFE in the vacuum case are characterized by two families of multipole- 
moments: mass- and spin-moments, see for example ns mi m [ 25 ]. Some of these solutions 
might be appropriate to represent stellar objects where the field moments can be related with 
corresponding body moments as integrals over the field generating sources. Moreover, the post- 
linear approximation of these metrics must be compatible with the HT solution. Quevedo et al. 
and Frutos et al. compared the HT solution with exact and approximate solutions of the EFE with 
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a quadrupole moment Q mmm of first order. Comparisons with the second order in Q of these 
solutions are still missing. 

Geroch and Hansen (GH) defined a procedure to find the field multipole moments of static and 
stationary spacetimes mm- Alternative definitions of relativistic multipole moments were given 
by Simon and Beig [28] and by Thorne [29]. It is important to mention that the GH multipole 
moments m m are equivalent to the Thorne moments for stationary systems [H|. Using the 
Ernst formalism EE], Fodor et al. found an elegant method to find the multipole moments of 
a given spacetime m- This method was later generalized by Hoenselaers and Perjes [19]. The 
relevance of taking the correct relativistic multipole moments of numerical spacetimes for modelling 
astrophysical objects such as neutron stars was discussed by Pappas and Apostolatos [23], who used 
a method due to Ryan m to derive the multipole moments. 

Nowadays, the use of harmonic coordinates is costumary, since the form of the transfomed metric 
tensor using harmonic coordinates has a special structure that allows to read off the Thorne mo¬ 
ments directly even for the non-stationary case [181129] , The multi-polar post-Minkowskian formal¬ 
ism (MPM) that was developed by Blanchet, Damour and Iyer [21 [3 is also formulated in harmonic 
coordinates. Applying this formalism Blanchet found spacetimes containing mass-quadrupole and 
quadrupole-quadrupole terms lam- The main goal of this paper is to compare these results with 
the ones we get from a static HT approximation with squared quadrupole moment. 

The paper is organized as follows. In the second section, we briefly describe the HT metric. We 
find a new expanded version of the HT metric with a squared quadrupole moment, in the third 
section. In the fourth section, the harmonic coordinates for this HT solution are obtained and the 
metric is expressed in these. In the fifth section, it is shown that our HT harmonic metric and 
the Blanchet metric coincide at our level of approximation. Finally, a coordinate transformation is 
found from the post-linear version of the Erez-Rosen metric [510 [301 EU [32] to this HT solution 
with no rotation in the sixth section. 


2 The Hartle-Thorne Metric 


The Hartle-Thorne metric [161 26] is an approximate solution of vacuum EFE that describes the 
exterior of any slowly and rigidly rotating, stationary and axially symmetric body. The metric is 
given with accuracy up to the second order terms in the body’s angular momentum, and first order 
in its quadrupole moment. It therefore has three parameters: mass M, spin S and quadrupole- 
mornent Q. The HT solution is given by 

ds 2 = g tt dt 2 + g rr dr 2 + g ge dQ 2 + gee sin 2 9dcj) 2 + g t(j) dtd(j), (1) 

with metric components 
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The functions Q 2 are associated Legendre polynomials of the second kind 
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3 The Post-linear Hartle-Thorne with Q 2 term 

Neglecting m 3 , J 2 -terms and changing q —>• — q in the HT-metric one obtains 
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is solution of the EFE by means of a REDUCE program m- At this point, it is possible to add 
rotation into this metric to first order without problems. 
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4 Transformation to Harmonic Coordinates 


Harmonic coordinates (cT, X, Y. Z) are especially useful, because the form of the metric tensor in 
these coordinates has a special structure where the Thorne-moments can be read off directly also 
for the non-stationary case [18, 29]. The harmonic coordinate condition reads 


The solution of this equation for the HT metric including g 2 -terms reads 
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The transformation of the metric from the (r, 9, 4>) coordinates to such harmonic coordinates 
(. R , (j)), was performed by means of the differentials 1-forms (dX, dY , dZ). From these 1-forms, 
it is solved for the other 1-forms (dr, d9, dcj >) 
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Substituting (11111 . (1121) . and (1131 into the metric with metric components (O, the metric in Cartesian 
coordinates to post-linear order takes the form 
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5 Comparison with known results 

Our metric (11511 can directly be compared with the one derived in Blanchet [5[ 0) that was derived 
within the MPM-formalism which works with a tensor field 
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and uses the Landau-Lifshitz form of the field equations in harmonic gauge (see e.g., 0] for more 
details). 

The post-linear metric components under this harmonic coordinates can be written in the following 
form 
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where Si is the total angular momentum of the object, and rii = x l /r with r 2 = x 2 + y 2 + z 2 . In 
our case, we only have the z component, i. e. S z = male. 

The mass-quadrupole and quadrupole-quadrupole metric components as obtained by Blanchet mm 
take the form 
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and are the symmetric and trace-free parts (e.g., [2H [29] ) of the Cartesian tensor 


ni 1 ...i l =n h ---n il . (22) 

For our axially symmetric body the quadrupole-moment Q appears in the Cartesian quadrupole 
mass-tensor in the form 


Mij = ~(3ij ~ 3<5i3<5j3). (23) 

Substituting these M ab components into (1211) and the resulting h a p into (1201) . we get the metric 
components lfT5jl . 

6 Transformation from Erez-Rosen to Hartle-Thorne 

The Erez-Rosen (ER) metric represents an static exact solution of EFE with axial symmetry and a 
quadrupole moment [5U8U30U3I1132] . Keeping only q 2 terms the ER-metric in spherical coordinates 
(ct, r, 6 , <))) reads 


9tt 


-[1-2U- ±^p 2 - ±,U*P 2 + JL q Wp? 


7 



9rr 


= l + 2U + AU 2 + 4 qU 3 P 2 + 4<7t/ 4 (5P 2 2 + 11 P 2 - 1) 

15 45 

+ ^^ 6 (25P 2 3 - 12 Pi - 6 P 2 + 2) (24) 

g ee = r 2 (l + ^qU z P 2 + ^qU 4 (5P 2 + 5P 2 - 1) 

+ ^9 2 t/ 6 (25P|-12P 2 2 -6P 2 + 2 )) 

g H = r 2 sm 2 e^l + ^qU 3 P 2 + ^qU 4 P 2 + ^q 2 U 6 P 2 y 

where 


Q ^ rr m 

q = —and 17 = —. 

m 6 r 

The following transformation converts the ER truncated metric into the static HT metric at the 
same level of approximation 


r 


6 


R 


mg g 

1 + lP fl + R^ h 


0 


mq q 2 

+ IP 91 + RS 92 


(25) 


where 


h = -L(10P 2 -8P 2 -2) 

/2 = i ^(40P|-24R|-43) 

gi = — — (5P 2 — 2) cos 0 sin 0 
45 

92 = ^ 5(12 “ 30P 2 ) cos 0 sin© 

with P 2 = P 2 ( cos0). 

The transformed metric components are given by 


9tt — 
9R& - 

9RR = 

+ 


gee = 


- I 1 - m - |^ 3p2 - §^ 4p2 + W ^ P2 

0 

1 + 2U + m 2 + 4 qU' i P 2 + \gu 4 p 2 
45 3 


1 


—q 2 U & ( —P 2 - —P 2 + 

ni r i l a r Z a rr ^ 


315 


168 


336 


45 


12 


45 

30 


1617\ 
~l5~ ) 


R A 1 + ir=qU a P2 + -rzqWPi 


45 


q 2 U 6 


315 


45 

308 o 56 „ 301 

^ + 45^ 16- 


(26) 


+ 



12 


30 


= R z sin" 0 1 + —qU 6 P 2 + —qU i P 2 


45 


45 


1 


315 


308 


56 


301 


+ U ( ~JT P 2 + 77^2 - -rzr 


45 


45 


45 


where U = m/R. 

Changing q —>• 15g/2, one obtains the static HT metric at this level of approximation. 


7 Conclusions 

We expanded the HT-metric and kept only linear terms in the rotation parameter and quadratic 
terms in the mass parameter. Then we included second order terms in the quadrupole parameter 
by solving the EFE in vacuum perturbatively. We then shows that this form of the metric agrees 
with a corresponding metric that was derived within the MPM-formalism by Blanchet at the same 
order of approximation. 

A transformation linking our static HT solution with the ER metric expanded in Taylor series was 
also found. This provides a validation of all these metrics. The quadrupole moment is an important 
feature which is included as a physical parameter in all these solutions. 

These spacetimes can be used to represent realistic solutions of self-gravitating (axially symmetric) 
mass distribution of perfect fluid. This is because the HT solution can be smoothly matched with 
interior perfect fluid solution with physically reasonable properties. 
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